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This document gives a technical description of a multi-channel parametric 
audio coding system as developed by Philips. The goal of this system is to describe an m- 
channel signal by an n-channel signal, with n<m, and parameters describing a spatial image 
in order to reconstruct the mchannel signal. Although the techniques described in this 
document could be extended to coding any m to any n channels, the embodiments described 
in this document is to provide a technical description of coding 5(.l) to 2 or 5(.l) to 1 channel 
coding. The extension ".1" denotes the presence of an LFE channel. Furthermore, it is 
assumed that when reproducing the multichannel signals, a typical loudspeaker setup is used 
consisting of a Left front (Lf), a Right front (Rf), a Centre (Cf), a Left surround (Ls), a Right 
surround (Rs) and optionally a low-frequency effects (LFE) speaker. 

HIGH LEVEL DESCRIPTION 

Figure 1 shows a general block diagram of the multi-channel parametric 
encoder. The multi-channel input signal consisting of the five channels Lf, Rf, Cf, Ls, Rs and 
the optional LFE channel are analyzed resulting in a set of parameters describing the spatial 
image. Depending on the configuration either a mono down- mix channel M is generated or a 
stereo down mix consisting of the left and right channels Ld and Rd is generated. This mono 
(M) or stereo (Ld, Rd) signal is then encoded using a conventional mono or stereo audio 
encoder respectively. The bit stream resulting from this encoding process is merged with a bit 
stream derived from the coded spatial parameters preferably in a backwards compatible 
fashion. 
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Figure 1 : Block diagram of generalized multichannel parametric encoder. 
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A block diagram of the corresponding multi- channel parametric decoder is 
given in Figure 2. First the bit stream is de- multiplexed resulting in a (backwards compatible) 
bit stream for the mono or stereo audio decoder and a spatial parameter bit-stream. The mono 
or stereo decoder then reconstructs the coded mono down- mix signal M' or the stereo down- 
mix signal (Ld',Rd') respectively. Concurrently, the spatial parameters are decoded. Finally 
the multi-channel signal is reconstructed by imposing the spatial parameters onto the down- 
mix channel(s). 
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Figure 2: Block diagram of generalized multi- channel parametric decoder. 



DETAILED TECHNICAL DESCRIPTION 
Time/frequency transform 

The multichannel analysis and reconstruction blocks require spatial analysis 
and synthesis to be performed on individual time/frequency tiles. Therefore, a time/frequency 
transform is required with the following prerequisites: 

The transform is preferably complex, to enable measurement and modification 
of (relative) phase values between input and output channels; 

The transform should be oversampled, to avoid abasing distortion which 
would result from time and frequency dependent changes in a critically- sampled system; 
The frequency resolution should be non-uniform according to the frequency resolution of the 
human auditory system; 

The time resolution is generally rather low, except in the case of transients. 

A generalized block diagram of a spatial encoder is shown in Figure 3. A 
multi-channel input signal is first transformed to the frequency domain. Subsequently, a 
downmix and spatial parameters are generated. The downmixed signals are subsequently 
transformed to the time domain. The decoder basically performs the inverse process. 
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Figure 3: Block diagram of generalized encoder processing stage. 

Currently, two time/frequency transforms are used which meet the 
prerequisites mentioned above. The first transform comprises time-domain segmentation 
f oUowed by FFTs. All input signals are segmented and transformed by means of the Discrete 
Fourier Transform (DFT): 



N/2 



Xy [k] = £*,[» + * • A] • A, M • exp - j 



n=0 



2%kn 
N 



with *,[n] the i th time domain input signal, h a [n] the analysis window, / the frame index, h 
10 the frame update in samples, N the DFT lengfli and X iJ [fc] the DFT with frequency index 
k . At the output of the multi-channel encoder/decoder the processed signals are 
transformed back to the time domain by means of an inverse DFT: 

{n-i f 2%kn M 

gF u[fc ].exp|y— |, 

with Y u [k] the (zero-padded) DFT representation, y fi/ [n] the time- domain segment 
15 corresponding to frame I and h s [n] the synthesis window. The resulting output signal y, [»] 
is obtained by overlap-add of the segments y u [n] . 

The second transform which is of particular interest for memory and 
computational complexity reasons is a complex- exponential modulated filter bank. 

In the following sections, it is assumed that individual time/frequency tiles of 
20 all input channels are available for processing at both the encoder and decoder side. 

System based on Stereo downmix 
Encoder 

The aim of this encoder is to represent a 5.1 -channel input signal as a 
25 backwards-compatible stereo signal (i.e., with a spatial representation which resembles the 
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5.1 channel reconstruction as g)od as possible), combined with spatial parameters that enable 
reconstruction of a 5.1-channel output that resembles the original 5.1 input signals from a 
perceptual point of view. The structure of this system is depicted in Fig. 4, 
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5 Figure 4: Block diagram of an encoder providing a stereo downmix. 



10 



15 



The encoder consists of three parallel stages which all convert a stereo input 
signal to a mono signal, combined with parameter extraction which represents the spatial 
cues between the respective input signals. Each of these three parallel blocks computes 
(assuming input signals Xi and and output signal Xm): 

The ratio of the powers of corresponding time/frequency tiles of the input 
signals (which will be denoted Tnterchannel Level Difference', or ILD), given by: 



ILD =101ogl0 



) 



The average phase difference (or the phase difference which maximizes the 
correlation between the input signals), which is referred to as 'Interchannel Phase 
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Difference', or IPD, given by: 

ipd = ^£x y mx;vcij. 

The coherence (ICC), which is the normalized cross-correlation between the 
input signals given the IPD value as in (2): 



ICC = 



Each paraUel processing block generates a parameter set which comprises all 
or a selected number of these parameters for each time/frequency tile (depending on the 
desired parameter bitrate, some parameters are not transmitted). Besides parameters, each 
10 parallel stage also computes a single output signal, which is a linear (complex) combination 
of the two input signals: 
X m m = w l X l [k]+w 2 X 2 m, 

with wi and wz complex weights which depend on the extracted parameters (wi = f(IID, IPD, 
ICC)). Preferably each time/frequency tile of Xn has a power that is equal to the sum of the 
15 powers of the input signals Xi and X 2 . 

A fourth parameter which is required for reconstruction of phase differences in 
the encoder is the Overall Phase Difference (or OPD), which is defined as the average phase 
difference between the first input signal and the mono output signal. 

The next stage of the encoder performs a down- mix from three (L, C, R) to 
20 two down-mix channels (Lo, Ro), combined with corresponding spatial parameters. Each 
down-mix channel is a linear combination of the input signals L, R and C: 

L 0 [k] =a l Hfc]+a 2 J?[*] + a 3 C[fc], 
i? 0 [fc] = p\L[fc] + MCfc] + MW- 

The parameters a* and p\- are chosen such that a good stereo image of the 
stereo signal consisting of L 0 [fc] and *>[*] is obtained. One of the prerequisites for a good 

25 stereo image is that a 3 equals 63 - 

At the decoder, channels L, R and C are predicted using the two down- mix 

channels Lo and Ro as follows: 
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im = C UL L 0 lk] + C 2 , lt R 0 [k1, 
R[kl = C hR L 0 m + C 2yR R 0 [kl 
Cm = C lc L 0 [k] + C 2 , c R 0 \_kl 

To this end, parameters C 1>z and C 2 , z (for Z = L, R or C) are computed at the 
encoder and sent to the decoder. 

A minimal Euclidian norm of the difference of signal Z[k] and its estimation 
5 Z[k] is used as optimization criterion to find the parameters C 1)Z and C 2Z . The square of the 
Euclidian norm of the difference of the original input channel Z[k] and its estimation at the 
decoder Z[k] can be written as: 

k ' 

with 

10 Z[k] = C KZ L 0 [k] + C 2 , z R 0 [kl 

Minimization of £|z[£] -Z{kf leads to the following expressions: 
Ciz = < L 0 lk],Z[k] >* \\R 0 [kf- < R 0 [K],Zm ><L 0 [k],R 0 Uc] > 

IKwfl 2 Kw| 2 -|<^w^oW>| 2 ' 

C _ <RolklZ[k]> \\hlkj-<L 0 [klZlk]><L 0 [_k lR f) [k]> 
l L omflR 0 [kf-\<L 0 lklR 0 [k] 

with 



\\A[kf = £|A[*f , 

k 

< A[k], B[k] >= ^A[k]B*[k]. 

k 



15 For the coefficients C x ,z and C 2)Z the following relations can be derived: 
«iC 1 . i +oc 2 C ljR +a 3 C 1>c =1, 

«iC 2 , £ +a 2 C 2R +a 3 C 2c =0, 
^iC XL +^ 2 C^ + ^ 3 C uc =0. 

Having 6 variables (C hL ,C 2>L ,C 1>R , C 2;R ,d, 0 and C 2 , c ) and at the same time 4 
relations between these parameters, only 2 parameters need to be sent to the decoder. C hL 
20 (henceforth notated by P) and C 2 ,r (henceforth notated by y) are transmitted to the decoder 
because of their similar statistical distributions. 
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Before the actual downmix is applied (i.e., L, C and R are combined to Lo and 
Ro), the signals L, C and R are preferably pre-conditioned by first applying a phase shift to L, 
R, and/or C to assure that the signal pair L and C on the one hand and R and C on the other 
hand have a non negative correlation to minimize energy loss by the downmix process. The 
applied phase shifts can be transmitted without any bit-rate costs by superimposing these 
phase shifts on the OPD values resulting from the individual two-to-one stages. 

If no LFE channel is present at the encoder input, the LFE channel can be 
considered as containing zeros only and all related processing steps can be discarded. In that 
parameter set 2 as shown in Figure 4 is irrelevant and does not have to be transmitted. 



case, 



Decoder 

The decoder basically performs the reverse process as depicted in Figure 4. 
In a first stage, the stereo input signal (Lo, Ro) is converted to a three-channel signal (L, C, 
R) based on parameter set 4. The upmix, based on the transmitted parameters 6 and y, is 
performed as follows: 

L[fc]=pZo[fc]+(Y-l)tfoOT 

#[*] = (p-l)L 0 [fc]+Y*o[*] 
cm = (1- P)L 0 [fc] + (1 -y)i?„W 

If a 3-channel reconstruction is desired, the decoding process is finished. 

For e.g. a 3.1, 5 or 5.1-channel reconstruction, the mono signals L, C, and R 
are subdivided into two-channel signals, based on the spatial parameters corresponding to 
each signal. The general structure of the mono-to-stereo upmix block is given in Figure 5. 
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Figure 5: Generalized structure of the mono-to- stereo upmix stage. 



A local copy of the mono input signal X is fed through a decollation filter H. 
This filter can be implemented as a (frequency- dependent) delay or reverberation module. 
The mono input signal X and its decollated copy Q are subsequently combined in a mixing 
stage to form the stereo output signal (Ty, Y 2 ): 
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cos(T) 0 Tzm" 
. 0 sin(T)JL<2M 



x = arctan 



a = —arctan 
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//£> /20 



S = 10 

If desired, this process is repeated for all three signals L, C, and R, resulting in 
a 5.1 (Lf, Rf, Ls, Rs, Cf, LFE) signal. It should be noted that the decorrelation filter H can be 
different for different mono -to- stereo processing blocks. 



Enhancement methods 

The backwards-compatible stereo downmix will usually result in a 
significantly-reduced spatial image compared to the original 5-channel reconstruction. It is 
possible to enhance the stereo downmix in such a way that the perceived spatial image of the 
downmix resembles the 5-channel reconstruction more closely by introducing virtual 
surround loudspeakers (or sometimes denoted as 'stereo widening' algorithm). The 
generation of virtual surround loudspeakers is based on cross-talk cancellation principles. 
Various methods for stereo- widening are currently available. However, these systems have 
an important drawback: if they are applied on a stereo downmix they result in a widening of 
the complete sound stage instead of a widening of the surround signals only. To overcome 
this drawback, we propose a cross- talk cancellation algorithm mat is applied on the stereo 
downmix, in which the amount of cross-talk cancellation depends on the spatial encoding 
parameters. Using this method, (1) only signal parts that would have been reproduced by 
surround loudspeakers in a 5-channel setup are processed, and (2) the cross- talk cancellation 
algorithm can be inverted, which is a requirement for high-quality 5-channel reconstruction. 
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System based on mono-downmix 

Approach 1 

Encoder 

A straightforward approach for a 5.1-to-l encoder is depicted in Figure 6. The 
encoder consists of subsequent stereo-to-mono analysis blocks, which are identical to those 
used in Figure 4. 
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Figure 6: Structure of 5.1- to-one encoder. 
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10 This approach reduces the number of input channels pair- wise using stereo-to- 

mono blocks. Each block generates a mono signal and spatial parameters, of which a 
selection is transmitted. A recommended parameter selection includes: 



Parameter set 1: 
15 IID, ICC for each time/frequency tile; 

IPD, OPD for time/frequency tiles up to about 2 kHz. 



20 



Parameter set 2: 

IID, ICC for each time/frequency tile; 

IPD, OPD for time/frequency tiles up to about 2 kHz. 
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Parameter set 3: 

HD only for time/frequency tiles up to about 150 Hz. 

Parameter set 4: 
5 HD, ICC for each time/frequency tile; 

Optionally: IPD, OPD for time/frequency tiles up to about 2 kHz. 

Parameter set 5: 

IID for each time/frequency tile. 



10 



Decoder 



A corresponding decoder consists of the reverse process as depicted in Figure 6. The 
individual building blocks are equal to those described in detail in Section 0. 

15 Approach 2 

Instead of performing a spatial analysis and downmix on a pair- wise basis, this approach 
performs the downmix in a single stage based on multidimensional Principal Component 
Analysis (PCA). Assuming 5 complex- valued input signals X|, the complex- valued 
covariance matrix Sy of the incoming signals is given by: 

20 a-, 

k 

Note that Sy = S/ . In principle the non-diagonal elements of the covariance matrix (i.e., i not 
equal to j) are complex. However, it is possible to extract and transmit the complex angles of 
Sy separately, resulting in a covariance matrix containing elements equal or larger than zero 
only. We will describe a real- valued covariance matrix from this point on. 

25 11 is assumed that the input signals Xj consist of a rotation (R) of a set of 

orthogonal signals Y; (i.e., R 1 = R T ) 
X = RY . 

Given the fact that the individual signals of Y are orthogonal, the covariance 
matrix of Y, S y , is diagonal. Consequently, the covariance matrix of X, S x , can be written as: 
30 S x =RS y R T . 

Given the diagonal matrix S„ the above expression equals the 
eigenvalue/eigenvector decomposition of S x . This means that the rotation matrix R and the 
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eigenvalues (which are equal to the energies of the orthogonal signals of Y) can be obtained 
by an eigenvalue/eigenvector decomposition of the co-variance matrix S x . 

The general encoder approach then consists of: 
computing and applying (relative) phase parameters of the input signals in order to result in a 
covariance matrix which contains smaller imaginary parts; a possible technique is applying 
complex principal component analysis; 

computing the eigenvalue/eigenvector decomposition of the complex covariance matrix, 

resulting from the original or modified input signals; 

(inverse) rotation of the input signals to obtain the principal components; 

transmission of the first principal component only (i.e., the component corresponding to the 

largest eigenvalue) as mono output signal 

transmission of: 

either the rotation matrix R (e.g. in terms of its angular components) and the (relative) 
powers of each principal component, or 

the (real- or complex valued) covariance matrix itself in terms of IPDs, correlations and 
power ratios; 

the relative complex phase angles; 

preferably an overall phase shift of the principal component signal. 

The decoder consists of the following steps: 
computing the rotation matrix R (either from the direct transmission or an eigenvalue 
decomposition of the transmitted covariance matrix; 

computing of a set of orthogonal signals by means of decorrelating the mono input signal; 
scaling the orthogonal signals by means of the transmitted (relative) signal powers; 
generation of a five-channel output by rotating the orthogonal signals; 
re-instating the complex phase relationships based on transmitted (relative, IPD and overall, 
OPD) phase information. 

Approach 3 

The five-dimensional PCA as described as Approach 2 is relatively complex compared to the 
cascaded approach described in Approach 1 especially since it is not possible to derive the 
eigenvectors analytically. As a compromise between both approaches the cascaded approach 
could be combined with three-dimensional PCA. This is illustrated in Figure 7. 
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Figure 7: Hybrid encoder approach. 



The front (F), surround (S) and centre (C) channels are derived similar to 
Approach 1. Consequently these signals are used as input to a three-dimensional PCA 
procedure resulting in a 3x3 rotation matrix R. The same procedure for encoding and 
decoding as described in Approach 2 can be used. 

Alternatively, the three channels input to the three-dimensional PCA 
procedure could consist of 
Downmix of Lf and Ls 
Downmix of Rf and Rs 
Downmix of Cf and LFE 



Residual-coding extensions 

Most encoder processing blocks reduce the number of input channels, 
combined with a parameterization of the relations between the original input channels. In 
principle, each reduction in number of channels results in an output signal, but also in a 
complementary residual signal (or more residual signals), which would in principle allow 
perfect reconstruction of the original input signals. In the decoder, these residual signals are 
artificially regenerated by decorrelation filters. It should be noted, however, that it can occux 
for certain time/frequency tiles that this synthetic residual signal is inappropriate for a 
perceptually high-quality decoder output signal. For these time/frequency tiles, it is possible 
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to encode and transmit the actual residual signal from the encoder, while for remaining 
time/frequency tiles, the decoder can use the synthetic residual signal. The scalability options 
of this approach are of interest; the encoded residual signals can simply be removed from the 
bitstream. In that case, the decoder will default to its synthetic residual signal. This approach 
can be used for three-to-one, three-to-two, as well as two-to-one processing blocks. 



Bit-stream syntax 

As already mentioned on page 1, the basic techniques used for multi-channel 
parametric coding can be applied to code different mum-channel signal configurations. This 
10 is reflected in the proposed bit-stream syntax. Some example possibilities: 

encode an Lf, Rf, Ls, Rs (4 channels) multichannel signal into a mono or stereo-compatible 
signal, 

encode a L(f),C,R(f),LFE (3.1 channels) multi-channel signal into a mono or stereo- 
compatible signal or 

15 encode a Lf,C,Rf,Ls,Rs,LFE (5.1 channels) multichannel signal into a L(f), C, R(f) multi- 
channel compatible signal. 

The bit- stream should be very flexible for decoding subsets of the multi- 
channel signal. Consider the situation where the original signal consisted of a 5.1 muh> 
channel signal encoded into a stereo signal using a similar structure as depicted in Figure 4. 
20 Furthermore, assume that reconstruction is required for a 3. 1 setup, consisting of a Left 

(front), Right (front), Centre (front) and an LFE speaker. By only decoding using parameter 
sets 2 and 4 the 3.1 channel signal is already obtained. Hence, it is not necessary to further 
decode the surround channels. 

In the bit-stream syntax the flexibility described above is obtained by 
25 explicitly defining the channel configuration of each encoder/decoder step in the 

mc_channel_config element. In this element it is described which parameters belong to 
which (intermediate) input and (intermediate) output channels (see Table 16). 

Consider the very simple case where the input signal consists of the mono 
signal M and the output signal consists of the mono signal M' and the signal LFE. If the 
30 frequency range of the LFE signal is only limited wilh respect to the bandwidth of the mono 
signal M, which is typically the case, only for the lower part of the frequency range the 
signals M' and LFE are obtained. For the higher part of the frequency range the signal M' is 
simply obtained as M'=M. A similar situation might occur in a 1 to 3 decoding step. Consider 
the case where the input signal consists of mono signal M and the output signal consists of 
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the left signal L(f), the right signal R(f) and the LEE signal. For the frequency range for 
which the LFE has been included, parameters for all three signals are included, which may 
e.g. consist of ED values between L and R and between L and LFE and coherence values 
between L and R, L and LFE and R and LFE (3 sets!). For the higher frequency range, the 
parameters will then only consist of IID values between L and R and coherence values 
between L and R. For this frequency range the decoder operates in a 1 channel to 2 channels 
decoding mode. 

For the elements datalto3() two alternatives are given, the first one uses a 
representation containing the parameters of the covariance matrix, split into power ratios and 
coherence values. The second alternative uses an angular representation (e.g. as defined by 
Euler) of the rotation matrix R. 



Table 1 - Syntax of mc_data() 



Syntax 


Num. bits 


Mnemo 
nic 


mc_data(){ 

if (mc_channe!_config==l) { 

mc_channel_config() 

for (s=0; s<nr_steps; s++) { 
switch(method[s]) { 
case 0: 

datalto2() 
break; 

case 1: 

data2to3() 
break; 

case 2: 

dataltoSQ 
break; 

case 3: 

dataltoSQ 
break; 

case 4: 

datalto51() 
break; 

} 

} 

} 


1 


uimsbf 

Note 1 


Mote 1: nr_steps is denned m the mc_channel_config() element Hence, fc 
first instance of mc data() mc channel config should be set to %1 


>r the 



15 
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Table 2 - Syntax of mc_channel_config() 





Num 
.bits 


Mnemo 
nic 


mc_channel_config() { 

nr_downmix_ch_coded 
nr_output_ch_coded 


1 
3 


Note 1 
Note 2 


nr__steps = 0; 

nr Qcc_jcxi — nr__ciu w iiiijlla l^h, 






while (nr_dec_ch < nr_output_ch) { 

for (ch=0; ch<method_in_ch[nr_steps]; 

input_ch[nr_steps,ch] 

for (ch=0; ch<method_out_ch[nr_steps]; 


3 


Note 3 






ch++) { 

output_ch[nr_steps,chJ 

} 

method[nr_steps] 


3 


Note 3 


3 




nr_dec_ch += 
nr_increase_ch[method[nr_steps]] ; 






nr steps += 1 

} 






num_env_default 

\ 


2 


Table 15 


Note 1: nr downrnix_ch = nr_downmix_ch_coded + 1 

Note 2: nr_output_ch = nr_downmix_ch + nr_output_ch_coded; 

Note 3: See Table 16. LFE can never be input channel! 
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Table 3 - Syntax of datalto2 



Syntax 



datalto2{ 



Num. 
bits 



if (datalto2_header) { 
if (enable Jid) { 
iidjtnode 

} 

if (enable _icc) { 
icc_mode 

} 

if (enable_ipdopd) { 
ipdopdjnode 

} 

freq_res 

if (var_framing) { 

for (e=0; e< num_env; e++) 



env_pos[e] 



} 



for (ch=0; 
ch<method_out_ch[nr_steps]; ch++) { 

if 

(output_ch[nr„steps,ch]==6) { 



} 



} 



if (enable_jid) { 

for (e=0; e<num_env; e-H-) { 
iid_dt[e] 

iid__data(e,nr_bands_coded) 

} 

} 

if (enablejcc) { 

for (e=0; e<num_env; e++) { 
icc_dt[e] 

icc_data(e,nr_bands_coded) 

} 

} 



1 
1 

99 
• • 

1 

?? 
1 

99 



1 

2 



nr_bands_coded ?? 



Mnemonic 



uimsbf 



uimsbf 



uimsbf 



Table 17 

uimsbf 
uimsbf 

uimsbf 



Note 1 



uimsbf 



uimsbf 
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if (enablejpdopd) { 

for (e=0; e<num_env; e++) { 
ipd_dt[e] 

ipd_data(e,nr_ipdopd_coded) 
opd_dt[e] 

opd data(e,nrJpdopd_coded) 
} 

} 


1 
1 


uimsbf 
uimsbf 


Note 1: In case one of the output channels is the LFM channel only a pari ot the 
signal is coded, denoted with nr bands _coded. . 
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Table 4 - Syntax of dta2to3 



Syntax 



data2to3{ 



e++){ 



ch++) { 
{ 



if (data2to3Jieader) { 

if (enablejbeta) { 
betajmode 

} 

if (enable_gamma) { 
gamma_mode 

} 

if (enabIe_opd) { 
opd_mode 

} 

freq_res 

if (varjraming) { 

for (e=0; e<num_env_Jcr; 



} 



env_posJcr[e] 



} 



} else { 

num_env_lcr = num_env_frame 
} 

for (ch=0; ch<method_out_ch[nr_steps]; 
if (output_ch[nr_steps,ch]==6) 
nr_bands_coded 

} 

} 



if (enablejbeta) { 

for (e=0; e<num_env; e++ ) { 
beta_dt[e] 

beta_data(e,nr_bands coded) 

} 

} 

if (enable_gamma) { 

for (e=0; e<num_env; e++) { 
gamma_dt[e] 

gamma_data(e 5 nr_bands_coded) 
} 

} 



Num Mnemo 
. bits nic 



1 | uimsbf 

99 

• • 

1 | uimsbf 

99 



1 

99 



1 

2 



99 



uimsbf 



uimsbf 
uimsbf 

uimsbf 



Note 1 



uimsbf 



uimsbf 
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if (enable_opd) { 

for (e=0; e<num_env; e++) { 
opd_dt[e] 

opd data(e ,nr_opd_co ded) 

} 

} 


1 


uimsbf 


Note 1: In case one of the output channels is the LFB channel only a part ot the 
signal is coded, denoted with nr bands^coded. . . 
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Table 5 - Syntax of datalto3 (Alternative A) 



Syntax 



Num. I Mnemonic 
bits 



datalto3{ 



if (dataltoSJieader) { 
if (enable Jid) { 
iid_jnode 

} 

if (enable_icc) { 
iccjnode 

} 

if (enable_ipdopd) { 
ipdopd_mode 

} 

freq_res 

if (var_framing) { 

for (e=0; e< num__env; e++) 



env_pos[e] 



} 



for (ch=0; 
ch<method_out_ch[nr_steps]; ch++) { 

if 

(output__ch[nr_steps ? ch]==6) { 

nrJbands__coded__set2 



} 



} 



} 



if (enablejud) { 

for (e=0; e<num_env; e++) { 
iid_dt[e] 

iid__data(e,nr_bands_coded) 
iid_dt[e] 

iid„data(e,n.r_bands_coded_set2) 
} 

} 

if (enablejcc) { 

for (e=0; e<num_env; e++) { 
icc_dt[e] 

icc_data(e 5 nr_bands__coded) 
icc__dt[e] 



1 | uimsbf 

99 

• • 

1 | uimsbf 

99 



99 



i 

2 



1 
1 



1 
1 



uimsbf 



Table 17 

uimsbf 
uimsbf 

uimsbf 



?? Note 1 



uimsbf 
uimsbf 



uimsbf 
uimsbf 
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icc_data(e,nr_bands_coded_set2) 
icc_dt[e] 

icc_data(e,nr_bands_coded_set2) 
} 

} 

if (enable_ipdopd) { 
nr = 

min(nr_bands_coded_set2,nr_ipdopd_coded); 

for (e=0; e<num_env; e++) { 
ipd_dt[e] 

ipd_data(e,nr_ipdopd_coded) 

ipd_dt[e] 

ipd_data(e,nr) 
opd_dt[e] 
opd_data(e,nr) 



} 



} 



1 
1 
1 



uimsbf 



uimsbf 
uimsbf 
uimsbf 



Note 1: In case one of the output channels is the LFE channel only a part of the 
signal is decoded to three channels, denoted with nr_bands_coded_set2. The rest 
of the signal is decoded to two channels, i.e., nr_bands_coded_set2 is set to 
nr_bands_coded. Note that for the bands decoded to three channels there exxst two 
sets of nDs, three sets of ICCs, two sets of IPDs and one set of OPDs. 
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Table 6 - Syntax of data 1 to 3 (Alternative B) 



Syntax 



Num. 
bits 



Mnemonic 



datalto3{ 



if (datalto3_header) { 

if (enable__angles) { 
anglejmode 

} 

if (enablejpdopd) { 
ipdopd_jnode 

} 

fireq_res 

if (varjxaming) { 

for (e=0; e< num_env; e++) 



env_pos[e] 



} 



for (ch=0; 
ch<method_out_ch[nr_steps]; ch++) { 

if 

(output__ch[nr_steps ? ch]==6) { 

nrJbands_coded_set2 



} 



} 



} 



if (enable_angles) { 

for (e=0; e<num_env; e++) { 
angle_dt[e] 

angle_data(e,nrjbands_coded) 
angle_dt[e] 

angle_data(e 5 nr_bands_coded_set2) 
} 

} 

if (enablejpdopd) { 



1 
1 

99 
• • 

1 

99 



1 

2 



99 



uimsbf 



uimsbf 



Table 17 

uimsbf 
uimsbf 

uimsbf 



Note 2 



uimsbf 



uimsbf 
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nr = 

min(nr_b ands_co ded_set2 ,nr_ip dop d_co ded) ; 

for (e=0; e<num_env; e++) { 
ipd_dt[e] 

ipd_data(e,nr_ipdopd_coded) 

ipd__dt[e] 

ipd_data(e,nr) 

opd_dt[e] 

opd„data(e 5 nr) 



} 



1 
1 
1 



uimsbf 
uimsbf 
uimsbf 



Note 1: In case one of the output channels is the LFE channel only a part of the 
signal is decoded to three channels, denoted with nr_bands_coded_set2. The rest 
of the signal is decoded to two channels, i.e., nr_bands_coded_set2 is set to 
nr_bands__coded. Note that for the bands decoded to three channels there exist two 
sets of DDs, three sets of ICCs, two sets of IPDs and one set of OPDs. 



The datalto5() and datalto51() elements are straightforward extensions of the 
datalto3() element with an increased amount of parameter sets. 



Table 7 - Syntax of iid_data() 



Syntax 


Num. bits 


Mnemonic 


iid data(e, nr_iid_par) { 

if(iid_dt[e]){ 

for (b=0 ; b<nr_iid_par; b-H-) { 
iid_par_dt[e][b] = 
sa huff dec(huff iid dt[iid quant], bs„codeword); 

} 


*>9 99 


bslbf 


} 

else { 

for (b=0 ; b<nr Jid^par; b++) { 
iid_par_df[e][b] = 
sa huff dec(huff iid df[iid quant], bs_codeword); 

} 

1 

1 


99 99 


bslbf 
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Table 8 - Syntax of icc_data() 



Syntax 


Num. bits 


Mnemonic 


icc_data(e 5 nr_icc_par) { 

if (icc„dt[e]) { 

for (b=0 ; b<nr_icc_par; b++) { 
icc_par_dt[e][b] = 
sajiuff_dec(huff ice dt ? bs codeword); 

} 


99 99 


bslbf 




} 

else { 

for (b=0 ; b<nr_icc_par; b++) { 
icc_par_df[e] [b] = 
sajiuff_dec(huff ice df ? bs codeword); 

} 

} 

} 


99 99 


bslbf 




Table 9 - Syntax of ipd__data() 


Syntax 


Num. bits 


Mnemonic 


ipd__data(e, nr_ipd_par) { 

if (ipd„dt[e]) { 

for (b=0 ; b<nr_ipdopd_par; b++) 

ipd_par_dt[e][b] = 
sajiuff_dec(huff ipd dt,bs codeword); 

} 






99 99 
• • • • • • • 


bslbf 


} 

else { 

for (b=0 ; b<nr_ipdopd_par; b++) 

ipd_par_df[e][b] = 
sa_huff_dec(huff ipd df,bs codeword); 

} 

} 

} _J 






99 99 


bslbf 
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Table 10 - Syntax of opd_data() 



Syntax 


iNum. DltS 


1VJLI1CII1UI11C 


opd data(e, nr_opd_par) { 

if (opd_dt[e]) { 

for (b=0 ; b<nr_ipdopd_par; 

b++){ 

opd_par_dt[e][b] = 
sa huff dec(huff opd_dt,bs_codeword); 

} 


99 99 


bslbf 


} 

else { 

for (b=0 ; b<nr_ipdopd_par; 

b++){ 

opd_par_df[e][b] = 
sa huff dec(huff opd df,bs_codeword); 

} 

} 

1 


99 99 


bslbf 



Table 11 - Syntax of beta_data() 



Syntax 


Num. bits 


Mnemo 
nic 


beta_data(e, nr_beta__par) { 

if (beta_dt[e]) { 

for (b=0 ; b<nr_beta_par; b-hf) { 

beta__par_dt[e][b] = ! 
sa huff dec(huff beta dt,bs_codeword); 

} 

} 

else { 

for (b=0 ; b<nr_beta_par; b++) { 
beta_par_df[e][b] = 
sa huff dec(huff beta df,bs_codew>rd); 

} 

} 

) 


99 99 


bslbf 
bslbf 


99 99 





5 
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Table 12 - Syntax of gamma_data() 



Syntax 


Num. bits 


Mnemonic 


gamma_data(e, nr_gamma_par) { 
if (gamma_dt[e]) { 

for (b=0 ; b<nr_gamma par; 

{ 

gamma„par_dt[e] [b] = 
sa_huf f_dec (huff_gamma dt,bs codeword); 

} 


99 99 


bslbf 


} 

else { 

for (b=0 ; b<nr_gamrna par; 

b++){ 

gamma_par_df[e][b] = 
sa_huff_dec(liuff_gamma df 5 bs codeword); 

} 

} 

} 


99 99 


bslbf 


Table 13 - Syntax of angle_data() 


Syntax 


Num. bits 


Mnemonic 


angle_data(e, nr_angle_par) { 

for (a-0 ; a<nr_angles ; a++) 
if (angle_dt[a,e]) { 
for (b=0 ; 

b<nr_angle_par; b++) { 

angle_par_dt[a][e][b] = 
sa_huff_dec(huff_gamma dt 5 bs codeword); 

} 


99 99 


Note 1 
bslbf 


} 

else { 

for (b=0 ; 

b<nr_angle_par; b++) { 


99 99 


bslbf 


angle_par_df[a][e][b] = 
s a_huff_dec (huf f_gamma df,bs codeword); 

} 

} 

} 

} 




Note 1: nr_angles follows from used method! 
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Table 14: Dependencies of variable method 



method 


description 




mpfhoft ont ell 


nr_incrcase_eh 


0 


1 to 2 


i 


2 


1 


i 


2 to 3 


2 


3 


1 


2 


1 to 3 


1 


3 


2 


3 


1 to 5 


0 (fixed) 


0 (fixed) 


4 


4 


1 to 5.1 


0 (fixed) 


0 (fixed) 


5 


5 


reserved 








6 


reserved 








7 


reserved 









Table 15: num_env as a function of num_env_default 



hum env default 


num env 


o 1 


0 


1 


1 


2 


2 


3 


4 



Table 16: Channel description 



input_ch/ 
output ch 


description 


abbreviation 


0 


mono 


M 


1 


left (front) 


L(f) 


2 


right (front) 


R(f) 


3 


left surround 


Ls 


4 


right surround 


Rs 


5 


centre (or front) 


C(F) 


6 


low frequency effects 


UFE 


7 


surround 


S 



Table 17: nrbands and nrbands_coded as a function of freq_res 



freq_res 


nr bands 


nr bands coded (per default) 


0 


10 


10 


1 


20 


20 


2 


34 


34 


3 


reserved 
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Annex A Parametric multichannel a u dio coder with 2. 3. 4 and 5 channel playback 
compatibility 

This invention is a multi-channel extension of the basic principle described in 
WO03/090208-A1. The 5-channel content is downmixed into 2 channels combined with a 
small amount of parametric overhead which enables 5-channel reconstruction at the decoder 
side. Moreover, 2, 3, and 4-channel reproduction are also supported. 

The important features of the proposed coder are: 
transmission of two audio channels (which can be encoded using an arbitrary stereo audio 
codec) and are preferably obtained using principal component analysis on the left- front and 
left-rear pair on the one hand, and using a separate principal component analysis on the right- 
front and right-rear signal pair; 

- transmission of parametric overhead, which includes: 

- inter- channel level differences between left- front and left-rear channels; 

- inter-channel level differences between right- front and right-rear channels; 

- inter-channel coherence values between left- front and left-rear channels; 

- inter-channel coherence values between right-front and right-rear channels; 

- the power ratio between the centre channel and the sum of the powers of left-front, 
left-rear, right-front and right-rear channels 

Additionally, the inter-channel phase differences and overall phase differences 
between left-front and left-rear on the one hand, and right- front and right-rear on the other 
hand, may also be included in the parametric bit stream. 

The parameters described above are typically analyzed as a function of time 
and frequency (i.e., for a set of time/frequency tiles). 



Encoder 

Assume a five-channel audio signal Ifin], £•[«], rfin], r r [n], c[n\, which describe the discrete 
time-domain waveforms for the left- front, left-rear, right-front, right-rear and centre signals, 
respectively. These five signals are segmented using a common segmentation, preferably 
using overlapping analysis windows. Subsequently, each segment is converted to the 
frequency domain using a complex transform (e.g., FFT). However, complex filter-bank 
structures may also be appropriate to obtain time/frequency tiles. This process results in 
segmented, sub-band representations of the input signals (which will be denoted by Lf[k], 
L r [k], Rflk], R r [Ic], and C[k] with k denoting frequency index). 



PHNL04038 8EPQ 



29 02.04.2004 
As a first step, the relevant parameters between left- front and left-rear are 
estimated. These parameters include the level difference (JID L ), the (average) phase 
difference (IPD L ) and the coherence (ICCd'- 

^L f mL* f m^ 



IID L =10 log 10 



IPD L =Z 



k 



ICC L = 



Y,L f mL* r [k] 



10 



15 



This process is repeated for the right-front and right-rear channels, resulting in 

IID R , IPD r , and ICC R . 

The second step consists of a principal component analysis (PCA) of the two 
signals left-front (Lf) and left-rear (L r ). To be more specific, these two input signals are 
rotated in order to obtain a dominant and a residual signal (Qffl), using a rotation 
angle a which maximizes the energy of the dominant signal: 



J L f [k-\.ex$(K-OPD L )) j 



cosa sin a 

- sin a cosa J_Zv OT.exp( j(-OPD L + IPD L 



20 



Here, the angle OPD L denotes an overall phase rotation angle, while IPD L 
ensures maximum phase- alignment of the two signals Lf and L r . The rotation angle a can be 
derived from the IID L and ICC L parameters following 



1 t f 2gICC L 
a=-arctan -p—y- 

with 



g = 10 IIDL ' 20 

25 
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The signal Q[k] is subsequently discarded, and signal Y[k] is scaled by a scalar 
fi to obtain L{k~\ in such a way that L[k] has the same power as the power of Q[k] plus the 
power of Y[k] (i.e., the signal Q[k] is discarded while the loss in signal power is compensated 
for by scaling of Y[kJ). It can be shown that the required scale factor fi is equal to: 



with 



l+ViT 



(j. = 1+ 



4ICC L 2 -4 



This process is also repeated for the right-front and right-rear signal pair, 
resulting in signal R[fc]. 

The last step entails mixing the centre signal C[fc] in both L[k] and R[k], 
resulting in the stereo output signal pair L 0 uilk~\, Rourlk] : 







~L[k]+zcm~ 


_R OUT [k ]_ 




_R[k] + eC[k] 



Here, e denotes a weight that determines the strength of C[k] in the downmix 
(typically 0.707). A parameter IID C that describes the power of C with respect to the power 
of L and R is extracted: 



IID C =10 log 10 



£ 2 Y,Clk]C*[k] 



^wtfm+^Riw'm 



The process as described above is repeated for each time/frequency tile. 
Subsequently, the signals Lourik] and Rourtk] are transformed to the time domain and 
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combined with previous segments using overlap-add, resulting in the output signals IovtM 
and rouAn]. A schematic overview of the encoder is shown in Fig. Al. 



lr 



Segment 

And 
transform 



Segment 

And 
transform 



rf. 



Segment 

And 
transform 



Lf 



Lr 



PCA 
rotation 



Parameter to 
PCA angle 



Parameter 
analysis 



Rf 
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Setl 



R 



PCA 
rotation 



Parameter to 
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Parameter 
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^Parameter 
Set2 



Mixing 
And 
Parameter 
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Inverse 
Transform 
and OLA 
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Set3 



Fig. Al. Schematic overview of the encoder 

Decoder for stereo playback 

For stereo playback, the transmitted signals lovrYril and r 0 vAn\ are reproduced 
10 over the two playback channels, without further processing. 



15 



Decoder for 3-channel playback 

For 3-channel playback, the two received channels lourin] and r 0 ui\)i\ are 
segmented and transformed to the frequency domain. Subsequently, the output signals 
R[k] and C[k] are obtained as follows: 



w L L OUT 

W r R 0 ut 

} RC ±x OUT. 



W L C L OUT+ W ^ R 



with 

20 
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0.5 fcj 

k 

o 2 R =J^R[k]R*[k] 

k 

C 2+10~" Dc/l0 



Decoder for 5-channel playback 

For five-channel playback, first the 3-channel playback reconstruction as 
described above is performed, resulting in signals ZJfc], R[k\ and C[k]. The next step entails 
splitting L[k] in Lfi]c\ and L r [k], and sphtting R[k] in Rfik] and R r [k]. This splitting process is 
performed using the inverse PCA rotation as used in the encoder. The dominant Y[k] and 
residual Q[k] signal, which are required for the inverse PCA rotation, are obtained as follows: 



~Y[kf 




L[fc]cosy 






#WL[£]sin y 



Here, H[k] denotes an all-pass decorrelation filter to obtain a decorrelated 
version of L[k\. The angle ? is given by: 



f 



J = arctan 



Subsequently, the signals L£k] and LJk] are obtained using inverse PCA: 
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cos a 


_L r m_ 




sin a 



-anaTexp(jOPD L ) 0 
cosa 0 exp{jOPD L -IPD L )\ 



Y[k1 



This process is then repeated for the right channel. 



5 Decoder for 4-channel playback 

The decoder for 4-channel playback can be simply obtained by first decoding 5 channels (//, 
l r , c, r fi and r s ), followed by mixing of the centre channel (c) in front left and front right: 



lf,playback= If +0.101 C 

10 r f>p i ayb ack= r f + 0.707 c 

The factor 0.707 ensures that the total power of the centre channel is constant, 
independent of playback through the single centre speaker or as phantom source created by 
left front and right front. The surround channels remain unchanged (i.e., the signals are the 
15 same as for 5-channel playback). 
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Annex B Coding the low frequency eff e cts (LFE) channel in parametric multichannel audio 
coding 

The multichannel audio coder described in on the previous pages, does not 
support the coding of a low frequency effects (LEE) channel, which is incorporated in a 
standard commonly known as 5.1. 

The proposed coder incorporates the LFE channel. To this end additional 
parameters need to be sent to the decoder. Also, the same time, the 2-channel down- mix is 
modified to enable the decoding of the LFE channel. This is done in such a way that the good 
quality of the stereo image of the 2-channel down- mix is preserved. 

The parameters are typically analysed as a function of time and frequency (i.e., 
for a set of time/frequency tiles). The bandwidth of the LFE is typically limited to 
approximately 120 Hz. The proposed solution allows for a variable bandwidth of the LFE 
channel. 

A multichannel coder is described in Annex A. A schematic overview of its 
encoder is shown in Fig. 1. in the following, only the changes required for incorporating the 
LFE channel are explained. 

Encoder 

In Fig. B2, a schematic overview of the encoder incorporating the subwoofer 
channel is shown. 

Assume that c[n] and lfe{n\ describe the discrete time-domain waveforms for 
the centre and the LFE signal respectively. The signal Cs and the parameters from the block 
'parameter analysis' are obtained from signals c and Ife in the same way in which the signal L 
and the parameters from the block 'parameter analysis' are obtained from signals If and Ir, as 
described in Annex A. There is however a difference. Because of the low- frequency 
behaviour of the signal Ife, this is done only for a limited number of frequency subbands. To 
this end, a parameter describing the number of frequency subbands occupied by the signal Ife 
is incorporated in parameter set 4. In the remaining higher subbands, only the signal C is 
transmitted. The other parameters included in parameter set 4 are the level difference (ED) 
and possibly the phase difference (D?D) between the centre and the LFE channel. If the JPD 
is sent, also the OPD parameter needs to be transmitted. 
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Fig. Bl. Schematic overview of the encoder described in Annex A 



Decoder for 5 J -channel playback 

For 5.1-channel playback, the signals C[k] and LF£|XI are obtained from the 
reconstructed signal Cs[k] similar to the way in which signals Lf[k] and Lr[fc] are obtained 
fromL[fc] as described in Annex A. The values of the not transmitted parameter ICC are set 
to 1. If the IPD and OPD are not transmitted, they are set to 0. 
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Fig. B2. Schematic overview of the encoder including an LFE channel 
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Annex C, Reducins the inter-chann el interference 

The multi-channel audio coder described in Annexes A and B aims at: 
Approximating the multi-channel audio with two channels and parametric overhead. This is 

done for low bit rate reasons. 

Backwards compatibility with 2, 3, 4 and 5-channel reproduction systems. To 
this end a good quality of the stereo image of the 2 transmitted channels is required. 

The multi-channel audio coder described in Annexes A and B exhibits a 
significant amount of inter-channel interference, mainly due to the demand for a good stereo 
image of the 2-channel down-mix. The current invention significantly reduces the inter- 
channel interference. In this way, the quality of the reproduced multi-channel audio is 

significantly enhanced. 

The coder proposed in this Annex C represents N input channels by 2 down- 
mix channels and parametric overhead. In order to get the best possible reconstruction, in the 
sense of least- square-errors, of the N input channels at the decoder using only 2 channels, 
Principal Component Analysis (PCA) should be used. Perfect reconstruction of the N input 
channels at the decoder is only possible if all N channels from PCA are used. Drawback of 
PCA is the fact that no control can be exerted over the 2 down- mix channels. This means that 
the abovementioned requirement regarding Ihe good quality of the stereo image is not met 

when employing PCA. 

As in the case of PCA, also in the case of employing 2 down-mix channels 
that do have a good quality of the stereo image, a perfect reconstruction at the decoder is only 
possible when these 2 down-mix channels are extended with an appropriate set of N-2 
channels. As opposed to PCA, whose N channels are orthogonal so that the N-2 discarded 
channels cannot be predicted using the 2 down- mix channels, now the N-2 channels can - to 
some extent - be predicted from the 2 down- mix channels. The proposed coder exploits this 
predictability at the decoder. In order to do so, parameters need to be sent to the decoder. 
These parameters are typically analysed as a function of time and frequency (i.e., for a set of 

time/frequency tiles). 

As compared to Annex A or Annex B tte following changes are made in the 

proposed coder: 

- at the encoder: Ihe parameters required at the decoder have to be computed. 

- in the bit- stream: the parameters required at the decoder have to be included. 

- at the decoder: a parameter dependent up -mix from 2 to 5 channels is performed (in 
Annex A, a fixed up- mix is performed). 
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Encoder 

Assume an iV-channel audio signal, where zi[n], z 2 [n], z N [n] describe the 
discrete time-domain waveforms of the N channels. These N signals are segmented using a 
common segmentation, preferably using overlapping analysis windows. Subsequently, each 
segment is converted to the frequency domain using a transform (e.g., FFT). However, filter- 
bank structures may also be appropriate to obtain time/frequency tiles. This process results in 
segmented, sub-band representations of the input signals, which will be denoted by Z x \k\ 
2i\k\ ... , Z N [k] with k denoting frequency index. 

From these N channels, 2 down- mix channels are created, being Lo[k] and 
Rolk], which are also segmented, sub-band representations. Each down- mix channel is a 
linear combination of the N input signals: 



i=l 



* 0 [*]=£p\-z ( .[fc]. 



(=1 



The parameters a, and 0,- can be set on the basis of a certain criterion. In ID 695741, this 
criterion is the good stereo image of the steieo signal consisting of Utk] and jfypfc]. 

Perfect reconstruction of the N input channels at the decoder is only possible 
when these 2 down- mix channels are extended with an appropriate set of N-2 channels. As 
opposed to PCA, whose N channels are orthogonal so that the N-2 discarded channels cannot 
be predicted using the 2 most relevant channels, now the N-2 channels can - to some extent - 
be predicted from the 2 down- mix channels. 

If the N-2 discarded channels are denoted by C 0 ,m, then these channels are 
predicted from the two down- mix channels by: 

C 0J [k] = C Ui L 0 [k] + C 2i R 0 [k]. 

For choosing parameters C u and C 2 . various optimisation criteria are possible. We choose as 
an optimisation criterion the minimal Euclidian norm of the difference of signal C 0 ,m and its 
estimation C 0 i [k] . Parameters C u and C 2> . need to be sent to the decoder. 

It can be shown that the parameters C u and C 2>j are related to the parameters 
that are obtained when minimising the Euclidian norm of the difference of the original input 
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channel Zm and its estimation at the decoder Z t [k] ■ A coder that uses these latter 

parameters is further described. 

The square of the Euclidian norm of the difference of the original input 

channel Zf[fc] and its estimation at the decoder Z f [fc] can be written as: 

5 £|z,^]-i,[^ a , 

k 

with 

zm = c uz L o m+c 2A R 0 iki 

Minimisation of £|z,.[fc] -Z,[kjf leads to the following expressions: 

10 

<L 0 m,z i [k]> yt 0 {kf-<R 0 m,z l m > <Lom, goW > 
Cl - Zi ' \hikin^[kf-\<L 0 m,R 0 m>\ 2 ' 

< Rpikizm > \\L 0 m\- < A>[fc],z,[fc] >*< L 0 m,R 0 iki > 
C2 - z ' ' \L 0 mf\Romt -\<L 0 vkwk-\>\' ' 

with. 

\Wf =Z\A[kf, 

k 

<A[fe],5[Jk]>=£A[fe]5'W. 

k 

15 

For the coefficients Ci,z and C 2 ,zu the following relations can be derived: 

N 

/=1 



20 
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Having N channels, with 2 parameters per channel (C hZi and C 2 , Z i for the ?-th 

channel), but at the same time 4 equations describing relations between these parameters, 2N- 

4 parameters need to be sent to the decoder. 

The process described above is repeated for each time/frequency tile. 

Subsequently, the signals Lo[k] and Ro[k] are transformed to the time domain and combined 

with previous segments using overlap-add, resulting in the output signals ib[«] and r 0 [«]. 

Summarising, the encoder sends 2 down- mix channels, k[n] and r 0 [«], and for 

each time/frequency tile 2AT-4 parameters, that describe how to retrieve the input channels 

from the 2 down- mix channels, to the decoder. 

Decoder 

At the decoder side, for each time/frequency tile, first the coefficients C h2i and 
C 2 ,zi are computed for all N channels, using the 2iV-4 coefficients that are transmitted and the 
4 equations describing relations between the coefficients. Then each input channel ZgR is 
approximated by Z, [&] , with 

Zm = C uz L 0 [k] + C 2 ^R 0 lk], 

where L 0 [k] and Rolk] are the received 2 down- mix channels. 

Incorporation of the coder in the multichannel coder described in Annex A or 

Annex B 

A schematic overview of the encoder of the multi-channel coder described in 
Annex A and Annex B is given in Fig. Al and Fig. B2 respectively. The coder described in 
this Annex C, can be used to replace the block called "Mixing And Parameter extraction", 
that has as inputs the channels L, R and Cs and as outputs the channels Lq andifo and 
parameter set 3. In order to get a good stereo image of the 2 down- mix channels Lo and R 0 , 
they are chosen as: 

R 0 lk]=R[k]+Cs[k]. 

Because of the three input channels (hence N = 3), only 2N-4, or 2 parameters 
need to be transmitted to the decoder. It is advantageous to transmit 2 parameters that have 
the same range (e.g. C ltL and C 2 ,r), so that the same quantisation can be applied to them. 

At the decoder side, for 3 or more channel playback, first all 6 parameters (C hL 
, C 2 ,l , Cut , C 2Ji Ci >Cj and C 2 ,cs) are computed using the 2 transmitted parameters and the 
relations between the 6 parameters. For example, if C UL and C 2>R are transmitted, then it 
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follows that C 2 , L = CSu-1, C W = Ci, L -l, C lt o, = 1 - C u and C 2)Cs - 1 - C 2 ,«. Hie output 



signals £[&], i?[£]and Cs[fc] are obtained as follows: 

C UL L 0 [k-] + C 2X R o m 
C^L,[k1 + C^ R R o m 
C IfC I 0 [*] + C a , c R 0 [fc]_ 



"to " 













5 



Playback of 4- or 5-channels is explained in Annex A. 
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Annex D Improved stereo codins 

Traditional coding schemes, like e.g. MPEG-1 Layer m (mp3, [1]) employ 
stereo coding tools to improve the coding efficiency. One of these coding tools is known as 
Mid/Side (M/S) stereo coding or Sum/Difference stereo coding [2]. Using M/S coding a 
stereo signal consisting of a left signal l[ri\ and a right signal r[n] is coded as a sum signal 
m[n] and a difference signal sin] 1 : 



m[n] = r[n] + l[n] 
s[n] = r[ri] - l[ri\ 



For (almost) identical signals /[«] and r{ri\ this gives a large coding gain as 
the corresponding difference signal s[n] is close to being zero, whereas the sum signal 
contains practically all the signal energy. Hence, in this situation the bit rate required for 
coding the sum and difference signals is close to the bit rate required for coding only a single 
channel. 

Alternatively the mid- side coding process can be described by means of a 
rotation matrix: 



r m[nT 
s[n] 



cos 



— sin — 



v 



"A f 

(TO) (K 
COS — 



l 4 




It is obvious that the left and right signals have been rotated over an angle of 
% I A . This is illustrated in Figure Dl. The sum signal can be interpreted as a projection of 
the left and right samples onto the line l=r, whereas the difference signal can be interpreted ; 
a projection of the left and right samples onto the line l=-r. 



1 Usually the sum and difference are calculated as m[ri] = c- (l[n]+ r[«.])and 

s[n\ = c ■ (l[nj - r{n\) . For explanatory reasons, the constant c has been discarded and the 

sign of s[n] has been inverted. 
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Figure Dl: Rotation of the left and right signals /[«] and r[n] over an angle of % I A . 

In order to obtain a minimum signal power in the residual signal (i.e., 
maximum coding gain) for a wide class of input signals, the rotation angle needs to be 
variable. An improved signal mapping, applied in a sub-band coding system using a variable 
rotation angle, is outlined in [3, 4]. The following signal mapping is applied: 



cos(oc) 
-sin (a) 



sin 



cos' 




where ra*[n] and s'{ri] represent the dominant and the residual signal respectively and the 
angle a is chosen to minimize the power of the signal s'{n) . Due to the unitary rotation 
operation, the power of m'[n] is then maximized. This process is illustrated in Figure D2. 
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r m' 
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Figure D2: Left and right signal mapping using a rotation along an angle a 

Using regular M/S coding (as represented by m and s in Figure Dl), the signal 
5 illustrated in Figure D2 would still result in a residual signal with considerable energy. As 

such, the coding gain obtained by M/S coding would be marginal. However, using the 

variable rotation angle a as illustrated above, a very small residual signal can be obtained. 

Obviously, this rotation technique works particularly well when the left signal is 

approximately a scaled version of the right signal. 
0 Both the M/S coding technique (i.e., rotation with a fixed angle) as well as the 

variable rotation technique described above are typically not applied to the broadband signal, 

but rather to signals (or frequency domain representations) representing only a smaller part of 

the full bandwidth of the audio signal, as e.g. described in [3, 4]. 

5 Although the rotation technique as described in [3,4] eliminates much of the 

disadvantages of M/S coding it is still sub optimal for signals having a strong phase or time 
offset. This is illustrated in Figure D3. 
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Figure D3: Rotation of left and right signals along an angle a . 

The oval- like structure indicates a time or phase delay between 1 and r. 
Because the rotation implies a real- valued projection, the residual signal will still contain a 
significant amount of energy (although usually less than using regular M/S coding). 

In order to further reduce the residual signal energy it is proposed to extend 
the current signal rotation by employing complex- valued phase rotations to the left and right 
signal components. From this point, it is assumed that the left and right signals are 
represented by their complex- valued frequency domain representations Z[fe] and r[fc] . One 
method to obtain such signal representation is as follows. First the left and right time domain 
signal segments are windowed: 

l q \n\ = l[_n + qHYh\n\ 
r g [n] = r[n + qH]-h[n] 

where q represents the frame index (q = 0,1,2,....), H represents the hop-size or update-size 
and n = 0..X-1 where L equals the length of window h[ri\ . These windowed segments are 
then transformed to the frequency domain by means of a Discrete Fourier Transform (DFT): 
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N-l 



»=0 



r[k] = J^r g in]-Gxp \-j 



«=0 



2%kn 

N 
2%kn \ 

N 



where iV represents the DFT length (N>L). Because of symmetry in the DFT only the first 
N/2 + 1 points are preserved. Furthermore, in order to obtain energy preservation, the first 
5 DFT points are scaled: 

/[0]=/[0]/2 
r[0] = r[0]/2 



10 



The signal model is extended in the following way: 



cos(oc) sin(a) Ye^' o Yffl N 
-sin(a) cos(a) I 0 Jr[&] 



15 



As can be observed from the equation above, the real valued rotation (using a 
variable rotation angle a is extended with a complex- valued phase modification matrix. The 
angle <p 2 is used to rninirnize the energy of the residual signal by (phase-) rotating the right 
signal. The common angle <p l can be used to maximize the continuation of the signal over 
frame boundaries. 

After signal mapping/modification the dominant and residual time domain 
signals m[n] and s[n] are obtained by first applying the inverse DFT on the frequency 
domain representations m[k] and s[k] : 



\ N J 
^ ( .2%kn\ 



N-l 

s g [k]=J^s[n]-ex.p 

n=0 



N 



where the dominant and residual frequency domain representations m[k] and s[k] have been 

zero-padded to length N . The time domain signals are then obtained by means of overlap- 
add: 



PHNL040388EPQ 



47 



02.04.2004 



m[n + qH] = m[n + qH] + 29t(m g [»] • h[ri\\ 
s[n +qm = sin + qH] + 23t{s g [n] • h[n]} 



Alternatively, complex-modulated filter banks could be employed to obtain a 
5 complex- valued frequency domain representation. 

As an example, the following synthetic signal is mapped by the three different 

methods described above: 



10 



/[«] = 0.5 cos(0.32n +0.4)+ 0.05 - z L [n] + 0.06z 2 [»] 
r[n] = 0.25 cos(0.32n +1.8)+ 0.03 • zjn] + 0.05z 3 [n] ' 

where zjn], z 2 [n] and z 3 [«] are independent white noise sequences with unit variance. Part 
of the signals l[ri\ and r[n] are shown in Figure D4. 



15 
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Figure D4: Left and right signals l[n] and r[n] respectively 

Figure D5, Figure D6 and Figure D7 show the results for the M/S transform, 
the signal rotation over an (optimal) angle a and the rotation over both an (optimal) angle a 
and phase rotation as proposed in this ID respectively. In this particular example the angles 
a , 9! and (p 2 are fixed. In a typical embodiment, these parameters are both time and 
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frequency dependent. In each figure, the top panel represents the dominant signal (m[n]); the 
bottom panel shows the residual signal (s[n]). 

The M/S mapping, as shown in Figure D5, clearly does not increase the coding 
efficiency for this particular situation. As a matter of fact, the residual signal energy, i.e., the 
energy of the signal s[n] , is higher than the energy of the input signal r[n]. 




850 900 950 1000 1050 1100 1150 1200 
n 



Figure D5: Resulting signals after M/S mapping 

The mapping by means of applying the (optimal) signal rotation a as 
illustrated in Figure D6 does also not help for this particular signal. Only a negligible 
reduction is obtained. 
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Figure D6: Resulting signals after rotation 

Finally, the results of the extension of the mapping as proposed in this ID are 
5 shown in Figure D7. Here a clear reduction of residual signal energy is observed. 
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Figure D7: Resulting signals after signal anaphase rotation 



A block diagram of an encoder according to the invention is given in Figure 
10 D8. The left and right frequency domain representations I and r are phase rotated to obtain 
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a maximum coherence (angle cp 2 ) and 

an optimal signal continuation over time (angle <p t ). 

Consequently, the phase-rotated left and right signals are rotated over an angle 
a for maximal reduction of the residual signal error as described above. The parameters a , 
(pi and q> 2 are quantized and coded into the bit stream. The dominant signal m and the 
residual signal s can be coded by two independent conventional mono audio coders (or of 
course one dual mono encoder). Additionally, certain parts of the time- frequency plane of the 
signal s, not perceptually contributing to the final output signal, can be discarded in the time- 
frequency (t/f) selector unit. The overall bit stream is formed by merging the bit stream 
corresponding to the dominant signal m, the residual signal s and the parameter bit stream. 
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Figure D8: Block diagram of proposed encoder 
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Figure D9 shows a block diagram of the decoder corresponding to the encoder 
of Figure D8. First the bit stream is de- multiplexed into separate bit streams for the dominant 
signal, the residual signal and the parameters. The bit streams for the dominant signal and the 
residual signal are decoded resulting in the signals rti 1 and s\ Then the inverse rotation (-a ) 
is applied to obtain preliminary left and right signal representations. Finally the left and right 
signals, V and r' respectively, are obtained by applying the inverse phase rotations (-<p t and 
-<P 2 ). 
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Figure D9: Block diagram of proposed decoder 
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The invention can also be advantageously applied in combination with a 
parametric stereo coding system [5]. A general block diagram of a parametric stereo coding 
scheme is given in Figure D10. It is fairly similar to the block diagram of the proposed 
encoder except for the fact that: 

- no residual signal is being transmitted and 

- the angle a is not transmitted, but instead an ED value and a coherence value ? are 
transmitted. 

The I1D value represents the Inter-channel Intensity Difference, denoting the 
(frequency and time variant) intensity differences between the left and right input channels. 
The coherence value denotes the coherence, i.e., the similarity, between the left and right 
input channels after phase synchronization. The angle a can be derived from the ED and 
coherence value. 
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Figure D10: Block diagram of parametric stereo encoder 

A corresponding decoder block diagram is shown in Figure D10. It 
corresponds to the block diagram of Figure D9 except for: 

- the residual signal s* is now estimated based on the dominant signal m 9 by means of a 
de- correlation process D and 

- the amount of coherence between the left and right output signals is determined by a 
scaling operation. 

The scaling operation basically describes the ratio between the dominant 

signal and the residual signal. 

m' 



CD 
-O 



X 

E 
i 

CD 

~o 

E 

CD 
CO 



Decoder 



Scaling 



Q" 1 



Signal 
rotations 



Phase 
rotations 



9/ 



Figure Dll: Block diagram of parametric stereo decoder 



PHNL0403 8 8EPQ 



52 02.04.2004 
Figure D12 shows a block diagram of the parametric stereo encoder enhanced 
with residual coding. With respect to Figure D10, the only difference resides in the coding of 
(part of) the residual signal s. Which part of the residual signal is coded by Coder 2, is 
determined by the time- frequency (t/f) selector unit. 



Phase 
rotations 



Signal 
rotations 



m 



Coder 



t/f 
selector 



Coder 



no, 9 («) 



Q 



x 

=5 

E 
E 

CD 
-♦— > 
O) 

CO 



E 
ca 

CD 

-i — > 
5 



Figure D12: Block diagram of enhanced parametric stereo encoder 

A block diagram of the parametric stereo decoder enhanced with residual 
coding is shown in Figure D13. The bit stream is first de- multiplexed into separate streams 
for the dominant signal, the residual signal and the stereo parameters. The dominant and 
residual signals are then decoded by Decoder 1 and Decoder 2, respectively. Those 
spectral/temporal parts of the residual signal which have been coded are signalled either: 
implicitly, by detecting "empty" areas in the time- frequency plane or 
explicitly, by means of bits in the parameter bit stream. 

This information is applied in the de- correlation unit D and the Combine unit 
to fill the empty time- frequency areas in the decoded residual signal with a synthetic residual 
signal. This synthetic residual signal is generated by using the decoded dominant signal m' 
and the de-correlation unit D. For all other time- frequency areas, the (transmitted) residual 
signal is applied to construct the signal s'. Note that for these areas, no scaling is applied. 
Hence, for these areas it can be advantageous to transmit the angle a in the encoder instead 
of the HD and coherence values as the bit rate required for the single parameter a is smaller 
than the bit rate required for the HD and coherence parameters. However, transmission of a 
instead of HD and coherence values makes the system non-backwards compatible to the 
regular PS system. The subsequent stages of the decoder operate in the same fashion as the 
conventional parametric stereo decoder. 
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Figure D13: Block diagram of enhanced parametric stereo decoder 

The criteria for deciding which time-frequency areas of the residual signal 
need to be coded are perceptually motivated. If (a time- frequency area of) the residual signal 
s does not contribute to the audio quality of the final decoder output signal, or if (a time- 
frequency area of) the de-correlated signal forms a perceptually valid representation of the 
(corresponding time- frequency area of the) residual signal, it is not necessary to code the 
residual signal. 

By coding different time-frequency aspects of the residual signal in the 
encoder and by multiplexing the corresponding data into a scalable bit stream, the enhanced 
parametric stereo codec can be extended to a bit-rate scalable codec. In a scalable system 
where the layers in the bit stream are dependent, the coded data corresponding to the 
perceptually most relevant time- frequency aspects should be placed in the base layer, and the 
less important data moved to refinement, or enhancement, layers. In this case the base layer 
would consist of the dominant bit stream, a first enhancement layer would consist of the 
stereo parameters and a second enhancement layer would consist of the residual bit stream. 

When layers are removed from the scalable bit stream, and information 
regarding the residual signal is thus lost, the enhanced parametric stereo decoder can combine 
the decoded residual signal reconstructed from the data in the remaining layers with the 
synthetic residual signal in the manner described above to form a meaningful residual signal. 
Furthermore, if a decoder is not equipped with a second waveform decoder (for the residual 
signal), e.g. due to complexity restrictions, the signal could still be decoded, although this 
would result in a lower quality level. 

Further bit-rate reductions can be obtained by discarding the values for cp t and 
cp 2 in the bit stream. In that case, the decoder reconstructs the output signals 1' and r' using a 
phase rotation of zero. This method effectively exploits the lack of sensitivity of the human 
auditory system to high-frequency (inter- aural) phase information. 
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Annex E Post -processing of the individual multi-channel signals in t he stereo down-mix 

Many of the multi-channel audio coders as described herein before generate a 
2-channel down-mix to be compatible with 2-channel reproduction systems. When this 
down-mix is created the spatial impression of the original multi-channel mix is lost. The 
current invention makes improvement of the spatial image of the down- mix possible after 
creation, based upon the parameters as determined in the multi-channel encoder. Also other 
post-processing techniques on the individual multi-channel contributions are made possible. 

The proposed method makes a reconstruction possible of the multichannel 
mix that is not affected by the post-processing. Also post-processing in the decoder is 
possible for stereo playback as a user- selectable, without the necessity to determine the multi- 
channel signal first. 

Without post-processing the down- mix is comparable with the standard ITU 
down-mix. The proposed method however improves the down- mix significantly. This is a 
very important issue, because it is very probable that the quality of the down- mix will be one 
of the selection criteria within MPEG. 

The proposed method is able to determine the contribution in the down- mix of 
the original channels in the multi- channel mix with the help of the determined parameters in 
the encoder. In this way post-processing can be applied to specific channels of the multi- 
channel mix (for example: stereo- widening of the rear channels), whilst the other channels 
are not affected. The post-processing does not affect the final multi-channel reconstruction. It 
can also be applied for an improved stereo playback without the necessity to determine the 
multi-channel mix first. 

This method differs from existing post-processing techniques in that it uses the 
knowledge of the original multi-channel mix (the determined parameters). 

Encoder 

Assume an N-channel audio signal, where z t [h], z 2 [n~\, ,z N [ri\ describe 

the discrete time-domain waveforms of the N channels. These N signals are segmented using 
a common segmentation, preferably using overlapping analysis windows. Subsequently, each 
segment is converted to the frequency domain using a complex transform (e.g., FFT). 
However, complex filter-bank structures may also be appropriate to obtain time/frequency 
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tiles. This process results in segmented, subband representations of the input signals (which 

will be denoted by ZJfc], Z 2 [£] Z N [lc\ with k denoting the frequency index). 

From these N channels, 2 down- mix channels are created, being L a [k] and 
R 0 [k] . Each down-mix channel is a linear combination of the N input signals: 

r=l 

The parameters a . and 0. are chosen such that the stereo signal consisting of 
L Q [k] and R Q [k] has a good stereo image. 

In the decoder the N input channels are reconstructed as follows: 

zm^c uz L 0 m^c 2tZ R 0 m 9 

where Z,[*] is an estimate of Z.[k] . The parameters C 1Zf and C 2 ^ are determined in the 
encoder and transmitted to the decoder. 




Figure El: The positioning of the post-processing and inverse post-processing block in the 
mult i- channel coder. 
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Figure E2: Basic scheme for post-processing of the stereo down- mix. 

On the resulting stereo signal post-processing can be applied in a way that it 
5 mainly affects the contribution of Z ( .[fc] in the stereo mix. In Figure El the position of this 

block in the codec is shown. 

Figure E2 shows how this post-processing block will look. The parameter w l 

determines the amount of post-processing of L Q [k] and w r of R 0 [k] . When w, is equal to 0, 

L 0 \Jc\ is unaffected, and when w t is equal to 1, L 0 {k] is maximaUy affected. The same 

10 holds for w r with respect to R Q [fc] . 

The following equations hold for the post-processing parameters w, and w r : 



15 



The blocks 2^ # 4 in Figure El are filters, which can be various types of 

filters, for example stereo widening filters (as shown in the end of this section). The resulting 
outputs are: 
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J^Ow J 



= H 



with: 



H = 



1- w l + \v l H l 
w,H 2 



l-w r + w r H 4 



When the filters H l H 4 are chosen properly, then the matrix H is 

invertible, when the filters are known in the decoder, because the parameters w, and w r can 
be calculated from the transmitted parameters. So the original stereo signal will be available 
again which is necessary for decoding of the multi-channel mix. 

Another possibility is to transmit the original stereo signal and apply the post- 
processing in the decoder to make improved stereo playback possible without the necessity to 
determine the multi-channel mix first. 

Incorporation of the coder in a multi-channel coder described in Annex A, B or C 

On these pages an encoder is described that codes 5-channel audio. The 
following equations are applied: 

L Q [k]=L[k] + Cs[k] 
R Q ik] = R[k] + Csm, 

in which Cs[k] is the mono signal that results after applying OCS between the LFE- 
(subwoofer) and center-channel. For L[k] and R[k] the following equations hold: 

Hk] = C, (cos(a,)L f +sin( a, )e jIPD 'L s ) 
R[k] = C r (cos(a r )R f + sin(a r )e jIPD 'R s ), 

where L f is the left/front, L s the left/surround, R f the right/front and R s the right/surround 
channel. The parameters IPD, and IPD r (inter-channel phase differences) and C, and 



10 
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C r (complex scaling parameters) are parameters that are determined in the OCS encoder. The 
angles in these equations can be calculated from the inter-channel intensity differences (IID) 
and the normalized cross-correlation (IC): 



2/C 10 //D,/20 
a, = 0.5 arctan( l miva _ x ) 

2IC lO IIDr ' 20 
a r = 0.5 arctan( /c Ia> r /io_ 1 )• 



In the decoder the following reconstruction is done: 

Lm = $L 0 [k}+<x-y)R 0 m 
kk] = (p -r)L 0 m+yRo[k} 

C[k] = (1- $)L o m + (1 -y)Ro [*]. 



where L[k] is an estimate of L[k] , R[k] an estimate of R[k] and C[k] an estimate of 
Cs[k] . The parameters (3 and y are determined in the encoder and transmitted to the 
decoder. 

15 Knowing all this, the functions that are used for the post-processing are: 

w t =/ 1 (a,)/ 2 (P) 
w r =/ 3 (a r )/ 4 (y), 

here f x f 4 can be any function. For example to apply stereo widening on the rear 

20 channels: 

/ 1 (a)=/ 3 (a)=sin(a) 

sin(0.57tp) if 0<|5<1 

/ 2 (P) = / 4 (P)= 1 if P^l 

0 if P<0. 



For the filter functions H l H 4 the following functions are then chosen (in the z-domain): 
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H x (z) = H 4 (z) = 0.8(1.0 + 0.2Z- 1 + 0.2 Z ~ 2 ) 
H 2 (z) = H 3 (z) = 0.8(-l .Oz' 1 - 0.2z~ 2 ). 

This invention can be applied in any multi-channel audio-coder that creates a 
compatible stereo down- mix. The invention can be applied in two ways: 
It can be applied before transmission to provide an improved down- mix for decoders that can 
only decode the stereo audio. 

It can also be applied in decoders that can handle the whole bit- stream as a 
user-selectable option to make improved stereo reproduction possible. 



PHNL0403 8 8EPQ 

6 1 02.04.2004 

CLAIMS: 



1_ An audio encoder as described hereinbefore. 

2. An audio decoder as described hereinbefore. 

53 An audio encoding method as described hereinbefore. 

4 An audio decoding method as described hereinbefore. 

5 An audio signal as produced by the method of Claim 3. 

10 

6 a storage medium having stored thereon a signal as claimed in Claim 
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ABSTRACT: 



This document gives a technical description of a multi-channel parametric 
audio coding system. The goal of this system is to describe an m-channel signal by an n- 
channel signal, with n<m, and parameters describing a spatial image in order to reconstruct 
the m-channel signal. 



Fig. 1 
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